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Abstract. Thermodynamic limit evolution of a closed quantum Heisenberg- 
type spin model with mean-field interactions is characterized by classifying all 
the symmetries of the equations of motion. It is shown that parameters of the 
model induce a structure in the Hilbert space by partitioning it into invariant 
subspaces, decoupled by the underlying Bogoliubov-Born-Green-Kirkwood-Yvon 
(BBGKY) hierarchy. All possible partitions are classified in terms of a 3 X 3 matrix 
of effective, thermodynamic limit coupling constants. It is found that there are 
either 1, 2, 4, or O(N) invariant subspaces. The BBGKY hierarchy decouples into 
the corresponding number of anti-Hermitian operators on each subspace. These 
findings imply that equilibration and the equilibrium in this model depend on the 
initial conditions. 



PACS numbers: 75.10.Jm,76.60.Es,05.20Dd 
1. Introduction 

Recent developments in ultracold atom and ion experiments have spurred renewed 
interest in dynamics of closed quantum many-body systems [TJ [5]. Although it 
is expected that equilibration (in a suitable sense) takes place in a generic non- 
integrable quantum system [3J, a substantial body of recent research has shown that 
the interim evolution, which is becoming accessible experimentally, is not yet fully 
understood [4]. In particular, several qualitatively different scenarios of equilibration 
that depend on the initial state, including relaxation to a state that is different from 
the thermodynamic equilibrium, have been reported in a closed non-integrable spin 
chain [S|. In this article a different spin model is studied, where the equilibrium state 
may not be independent of the initial conditions. 

A model, where N spin- 1/2 particles on a lattice interact by a mean- field Curie- 
Weiss potential, has been discussed in [5]. From a previous study of a related model 
[7], it had been known that, with a restricted choice of coupling constants and initial 
conditions, certain observables approach their microcanonical expectation values with 
a rate to that is system size-dependent and, in fact, diverges as t ~ N 1 / 2 . With 
the aim of extending the range of parameters, and better understanding the course 
of equilibration of homogeneous states with arbitrary initial conditions, a study of 
the general Curie- Weiss quantum Heisenberg model was undertaken in [6]. For that 
purpose, reduced density operators -Fi... n for all n = 1...N spins were expanded 
in terms of a set of tensorial coefficients {f n }n=o> anc ^ the evolution equations 
for all f n were derived, starting from the Bogoliubov-Born-Green-Kirkwood-Yvon 
(BBGKY) hierarchy. Remarkably, these evolution equations, expressed in scaled units 
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of time r = t/t have a non-trivial N — oo limit. This thermodynamic coefficient 
BBGKY hierarchy implies that in thermodynamic limit. i is the slowest time scale 
of equilibration for arbitrary initial conditions and parameters of the model. 

However, scaling is not the entire story, and additional information about 
equilibration is contained in the solutions of the equations which, unlike scaling, 
depend on a 3 x 3 matrix of spin-spin coupling constants J and on the magnetic 
field vector h in a non-trivial way. For example, it was shown that if the Curie- Weiss 
potential is of the form V = -NJ ± {S x S x + S y S y )-NJ z S z S z [ which corresponds to a 
diagonal matrix J = diag( Jj_, J±, J z )\, and where S a is the a £ {x, y, z} component of 
the average spin per particle, the coefficient BBGKY equations simplify considerably. 
Evolution of a certain subset of tensorial coefficients /„ could be determined 
analytically to all orders n = 1,2,..., showing a characteristic superexponential 
[~ cxp (— t 2 A/2), where A = (J z — J±) 2 ] approach to thermodynamic expectation 
values, which in this case are equal to zero. And yet, another subset of coefficients 
could be identified, which displays no relaxation at all. These analytic results are not 
expected to be qualitatively valid for more general couplings and, in fact, no such 
analytic results are known about as simple generalization, as J = diag( J x , J y , J z ), 
which differs from the one considered in [6J by allowing all three diagonal elements of 
J to be different. 

In this article, the thermodynamic coefficient BBGKY hierarchy equations are 
further investigated. The main goal is to classify, for a general 3x3 matrix J of 
coupling constants and a general magnetic field vector h = (h x ,h y ,h z ), all tensorial 
coefficients with different behavior. Although the number D(N) of such coefficients 
depends on the system size N and not on parameters J or h, it will be shown that 
the matrix of effective thermodynamic limit coupling constants J(J, h), which is a 
function of J and h, determines the partition of all the coefficients into several groups, 
decoupled by the BBGKY hierarchy. In this sense, J induces a topology in the vector 
space space of expansion coefficients. 

This study provides an insight that the -/V-body quantum dynamics may impose 
symmetries on the Hilbert space. It is believed in this work that similar symmetries 
are present in a class of more realistic lattice spin models, where interactions among 
particles depend on the distance between them. Existence of invariant subspaces 
implies that equilibration and the equilibrium state may not be independent of the 
initial state, but not only that. It allows to divide a large computational problem into 
several smaller ones and moreover, should be taken in consideration when approximate 
dynamical models, such as kinetic theories, are constructed. The latter include a 
closure scheme as one step of approximation, whereby an ad-hoc relationship among 
correlators, involving different numbers of particles, is postulated. The existence of 
invariant subspaces with rigorously independent evolution implies causal relations 
among spin-spin correlators which should be respected by a good closure scheme. 

The article is structured as follows. Section [5] contains the preliminaries: The 
Curie- Weiss anisotropic quantum Heisenberg model is presented, a special expansion 
of reduced density operators is recalled, and the derivation of the equations for the 
expansion coefficients is outlined [BJ. In section [31 a useful representation of the 
coefficient BBGKY hierarchy equations is derived. Moreover, since the subsequent 
analysis of the equations can be succinctly expressed using concepts of graph theory, 
part of section [3] is dedicated to the introduction to the relevant terminology. In 
section [4j equations are reformulated using the generalized Liouville operator. It is 
shown that this generalized Liouville operator is anti-Hermitian, provided that the 
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scalar product in the vector space is appropriately defined. In section [S] the partition 
of the coefficient vector space is investigated. It is shown that all the coefficients can be 
grouped in either four, two or 0{N) disjoint sets, depending on whether the matrix 
of effective coupling constants is diagonal isotropic or anisotropic, or has a certain 
number of non-zero off-diagonal elements. Finally, in section [5] conclusions are drawn. 

2. Model and a summary of previous results 

This section briefly outlines some of the previous results which will be necessary for 
the subsequent discussion. The Reader is kindly referred to [5] for details of derivation. 

2.1. Curie-Weiss anisotropic quantum Heisenberg model 

Consider a Hilbert space Km = <S>iLi c ? of N identical spin-1/2 particles, where the 
Hilbert space of each spin-1/2 particle is <Cf. The anisotropic Curie- Weiss Hamiltonian 
H\...n 1 acting on Hn, is defined by 

N N 
i—l i,j=l 

The local potential Hi and the interaction potential Vij are defined as 

Hi = - h *< . Vii = - X JabCT > b j ■ ( 2 ) 

aEX a,b£T 

where I = {x, y, z}, and where 

X = l/N. (3) 

The scaling factor A is necessary to render the energy per spin finite in the 
thermodynamic limit. The interaction potential Vij is determined by a symmetric 
3x3 matrix of spin-spin coupling constants J = (J a b). The local potential Hi is 
determined by the three components of the magnetic field vector h = (h a ). 

2.2. Operator BBGKY hierarchy 

The density operator of an TV-spin system is a self-adjoint, positive, trace-class 
operator S Hn- The normalization convention Tri. jv pn = 1 will be used, where 
Tri. jv denotes the trace over N degrees of freedom. The density operator satisfies 
the Von Neumann equation, associated to the Hamiltonian (JlJ 

ihdtPN = [Hx... N , p N ] . (4) 

Reduced n-particle density operators are defined by 

F\... n = Tr„ + i...iv pn 1 ^ n < N. (5) 

where Tr„ + i.. jv denotes a partial trace over N — n spin degrees of freedom. In this 
article, only homogeneous states will be considered. The pn for such states is invariant 
with respect to an arbitrary permutation of particle indices. As a result, the n-particle 
reduced density operators do not depend on which N — n particle complex is traced 
over in ([S]). 

Reduced density operators satisfy trace properties which are relations of the form 
Tr-n+i -Fl...n+1 = Fi...n- (6) 
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The set of evolution equations for all the operators {Fi n }^ =1 is called the 
Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy. It is derived from Q 
using the standard formalism [5]. For the Hamiltonian ([XJ) it reads 



J2 l V ihFi...n] + {N- n) Tr„ +1 ^ [V i , n+1 ,F 1 ... n+1 ] .(7) 



i<j 



The initial value problem for the reduced density operators i*i... n consists of the set 
of equations (J7) for each 1 ^ n ^ N, supplemented by the trace properties ([5]). 

It was shown [5] that there are two time scales involved in the dynamics: The 
slow time scale of long-range spin-spin interactions drives the equilibration and is the 
main focus of this article. This slow dynamics is coupled with the fast precession 
around the magnetic field vector h = (h x ,h v ,h z ), appearing in the terms Hi. The 
latter terms can be eliminated by passing to an interaction picture, defined by 



i=l 



Fi...„exp 



it 



Fi... n -> Fi...„ = cxp 
The BBGKY hierarchy in the interaction picture is then given by 

n n 
ihd t Fl... n = ]T A. 



i,3 = l 
i<j 



+ {N~n) Tr n+1 ^ [^,„+i(i), A-.n+i 
i=i 



where the two-body potential is now explicitly time-dependent, 

if 



Vij{t) = exp 









Vij exp 



- h {H i + Hj ) 



The interaction picture potential Vij (t) is of the same form as (J5|) , 

Vij(t) = -x J2 Mt)°i°i 



(8) 



(9) 



(10) 



(11) 



a.bel 



where the matrix J is replaced by a time-dependent J(t), defined by 

J(t) = B(th/h) T JB(th/h) . (12) 

Define b = th/h, h = h/ \h\, and c° = cos |26|. Moreover, define c u = h u sm\2b\, and 
c uv = 2h u h v sin 2 where u, v E {x, y, z}. Using these definitions, B(b) is written as 



+ c z 



B(b) 



c yy 



c° c 



(13) 



2.3. Coefficient expansion and scaled variables 

Convinced that a "good" coefficient expansion should automatically satisfy the trace 
properties for reduced density operators ([6]), and noting that the following definitions 



/ 1 °(o? + o5)+ ]T f 2 o-lo- 



F 



(14a) 
(146) 



a,bel 
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satisfy Tr 2 F\2 = F\ and Tri Fi = 1, in [S] we attempted a generalization by 

n 

F 1 ... n = 2-» j2Y,fs E ( 15 ) 

s=0aEX B pS<P 3 (n) 

where o~£ is the a-th component of the Pauli operator on the i-th lattice site. Here 
{/s}" = o i s the set of coefficients of expansion of F±,,, n , where f s £ j^3x...x3 j g a 
rank-s symmetric tensor with components /" , and where a = {a\02 ■ ■ ■ a s ) £ I s . By 
convention, fo = 1. Expansion (fTSj) holds also for in terms of coefficients {/ s }™ =0 - 
The set ^p s (n) consists of all s-element permutations (pi, ■ ■ ■ ,p s ) of particle 
labels pi £ {l,...,n} such that pi < pj for all 1 ^ i < j ^ n (i.e. all 
sequences are strictly increasing.) Considering for example s = 2, we have ^Pafa) = 
{(1, 2), (1, 3), . . . , (1, n), (2, 3), . . . , (n — 1, n)}. Lastly, a product of s Pauli operators 
on different lattice sites is denoted by 

s 

^ = li<- ( 16 ) 

i=l 

The expansion coefficients are related to the expectation values of many-particle 
spin observables by the identity 

fs = K> , (17) 
where (o* a ) is the expectation value of a product of s spin operators, defined by 

s 

<r a IV,.... (/ ,...,![< ). (18) 

i—l 

Mapping of the operator BBGKY hierarchy ([7]) and ^ onto equations for the 
coefficients {f s }^ =1 and {f s }^ = i is discussed in [5J. Here, only the final result for the 
interaction picture (|9]) is recalled: 

~dj n = A5-[/ n _i] + (1 - n\)v+[f n+1 ] , (19) 
where the components of are given by 

n 

[/„_!] = - E E ^bcJ6a i (*)/rr +C "°', (20a) 

b.cGX i,3=l 

n 

fe + ) (ai - a " ) [/„ + i] = - e Jut)J2 e ^J:+t i+c+d - (20b) 

6,c,d£Z i=l 

Here, e a f, c is the Levi-Civita symbol defined according to the convention e xyz = 1, and 

a — di + C—dj = (Ol, . . . ,Oi_i,C, Oj+l, . . . , Oj_l,Oj+l, . . . , On) (21) 

in (j20a[) is derived from a = (oi . . . a n ) by replacing the ith element by c and then 
deleting the jth entry of a, and 

a - a t + c + d= (oi, . . . , a,_i,c, o»+i, . . . , a„, d) (22) 

in (|20&p is obtained from a by replacing the ith element by c and then appending d. 

Tensors v^[f n ±i] are linear in the tensors f n ±\ and in the time-dependent matrix 
J(t). Moreover, v^; do not depend on A explicitly. A change of t and /„ variables to 

r = 2tX^ 2 /H, ~f' n = /n/A" /2 (21) 
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balances all the powers of A in the coefficients of in (|19p . These coefficients are then 
cancelled from the equation with the goal to arrive at a scale-free BBGKY hierarchy. 
The remaining dependence on the number of particles enters the scaled equations 
through the coefficient 1 — nX at order n, and in the matrix J(t). By (|TS|) (fT5|) . J(t) is 
a sum of a constant term J and time-dependent terms which, by the form of B(th/H), 
are trigonometric functions such as cos (2t \h\ /h). In r units of time, such terms are a 
fast oscillation on top of the mean J, and it was argued [6] that these fast oscillations 
can be taken into account by averaging (which becomes exact in thermodynamic limit.) 
With the additional assumption that /' varies slowly in the period of oscillation T\, 
f is replaced by its average value / (r) = (T A ) _1 f T /'(r')dr'. Application of this 
procedure on (|19p results in the scaled and averaged BBGKY hierarchy 

drf n « v- [/l_ x ] + (1 - An)T7+ [f' n+1 ] , (22) 

where are obtained from v„ by replacing J(t) by J. The latter is defined by 

J = T _1 / J(t)dt T = nh/\h\. (23) 
Jo 

The approximation in (|22j) is a result of a substitution /—>•/. It turns into an equality 
in two cases. Firstly, if there is no magnetic field in ([lj i.e. if h = 0, then there is no 
need to perform averaging. In this case the non-averaged coefficient BBGKY hierarchy 
has the same form as (|22[) where J is equal to the original coupling matrix J. Secondly, 
if separation of scales is a valid assumption^, the averaging is exact in the limit A = 0. 

This limit describes thermodynamic evolution of the model by means of an infinite 
dimensional set of differential equations. The structure of these equations is the main 
focus of this article. Dropping the bars and the primes from / and for simplicity, 
the thermodynamic limit BBGKY hierarchy is defined as 

d T f n = v-[f n - 1 ]+v+[f n+1 }. (24) 

The distinction between the original matrix J and the averaged one, J, is maintained. 

3. Desymmetrized coefficient BBGKY hierarchy 

This section introduces the desymmetrized coefficient BBGKY hierarchy and presents 
the framework for the subsequent study of its properties in sections U [5] 

3.1. Coefficient equations 

Equation (| 15[) maps F\...n onto a set of tensor coefficients {f n }n=oi where /„ is a rank- 
n three dimensional symmetric tensor with components and where a = (ai . . . a n ) 
is a coordinate multiindex, and a, G {x,y,z}. A symmetric rank-n tensor in three 
dimensions is specified completely by d(n) numbers, where 

n+l 

d(n) = £j = (n + l)(n + 2)/2. (25) 

3=1 



X It is justified rigorously in a special case, treated in [6] 
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The number of degrees of freedom D(N) is defined as the number of independent 
coefficients in the expansion of Fi jy. It is equal to a sum of the numbers of 
independent components from all the tensors /„ excluding the trivial f = 1, 

N 

D(N) = d{n) — (N + l)(N + 2){N + 3)/6 - 1. (26) 

n=l 

A sequence, formed from all the independent components of /„, can be viewed as a 
vector indep(/„) on a c?(rt)-dimensional vector space S n , 

indep(/ n ) G Sn ■ (27) 

where <S„ C R d ^ is called the ?vspin desymmetrized coefficient vector space, or simply 
the n-spin vector space. A sequence of all independent tensor components 

indep(/) = (indep(/i), . . . ,indep(/jv)) G Fn (28) 
can be viewed as a vector on a direct sum of n-spin desymmetrized vector spaces, 

N 

Fn = S n . (29) 

By the symmetry of f n , any /° is equal to /« , where 

a(m) = ( x, . , x , y, . . . , y, z, . . . , z) (30) 
m x elements 

is a permutation of a' with the components x, y, z arranged in contiguous blocks of 
m x ', m y , and m z labels. Hence, all independent components of /„ can be enumerated 
uniquely by a triple of nonnegative integers m S SDtn, where the set DJl n is defined by 

M n = {(m x ,m y ,m z ) :m a ^0,m x + m y + m z = n} . (31) 
In order to denote the allowed set of indices for an arbitrary n, define 

971 = {{m x , m v , m z ) : m a > 0,m x + m v + m z ^ N} . (32) 
Define the norm of m £ 9Jl„ by 



\m\ = m x + m v + m 



(33) 



As a result, all indep(/) G Fn are uniquely labeled by elements of 



N 



m (N) = \J Mn (34) 

n=l 

In the subsequent, instead of a cumbersome indep(/„), the sequence of all the 
independent components fn^ will be denoted by f n (m) G S n where m G 9Jt„. 

The coefficient hierarchy (f2~i|) is translated to a desymmetrized coefficient BBGKY 
hierarchy on Tn- This is formally achieved by replacing /° by f n (m) = fn , and 
by replacing v^[f n ±i] on the right-hand side of (|24[) by v^(m, f n ±i) 

d T f n (m,T) = u~(m,/„_i) + (1 - nA)w+(m,/„+i). (35) 

The derivation of v„ (m, / n ±i) is given in |Appendix A| It is convenient to separate 
the diagonal and off-diagonal components of J. Define vectors K = (K x , K y , K z ) and 
W= {W x ,W y ,W z ) as 

(K x i Ky , K z^) = \ J yy J zz i J zz «/ •, J xx «^ yy) j (36 fl) 

(W x ,Wy,W z ) = (j yz ,J zx ,J x y) . (36b) 
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Replacing f n ±\ by a generic / in (to, f n ±i) an d using the notation (|36a|) - (|366[) . 
terms u„ (to, /) for general K and W, and hence a general symmetric J, are given by 

v+(m, /) = - K x m x f{m x - 1, m» + 1, m z + 1) - K y m y f(m x + 1, m y - 1, m 2 + 1) 

- K z m z f(m x + l,m y + 1, to 2 - 1) 

- ^(m 9 - m^/fm 1 + l,™ 9 , m z ) - F^fra 2 - m^/fm 1 ,™ 9 + 1, m z ) 

- W^m* - m y )f(m x ,m y ,m z + 1) 

- m x W x [f(m x - 1, m y , to 2 + 2) - /(m* - 1, m v + 2, m 2 )] 

- m y W y [/(m* + 2, to 1 ' - 1, m z ) - f{m x , m v - 1, m z + 2)] 

- m z W t [f(m x ,m v + 2,m z - 1) - /(m x + 2, m a , m 2 - 1)] . (37a) 
v~(m; f) = K x m v m z f(m x + 1, m y - 1, m 2 — 1) + K y m z m x f(m x — 1, m v + l,m z - 1) 

+ K z m x m y f(m x - l,m y - l,m z + 1) 
+ W x m x (m y - m z )f(m x - 1, m y , m 2 ) + W y m y (m z - m x )f{m x , m y - 1, m 2 ) 

+ W z m z (m x -m y )f(m x ,m v ,m z - 1) 

+ raV 1 - ^[Wy/fro* - 2 im y + l,m z ) - W z f(m x -2,m y ,m z + 1)] 

+ m y (m y - l)\W z f{m x ,m v - 2, to 2 + 1) -W x f(m x + l,m y - 2,m 2 )] 

+ m 2 (m 2 - 1) [W x /(m x + 1, m y , m z - 2) - W y f{m x , m y + 1, m z - 2)] .(376) 

In (|37aj) - (|37b|) . indep(/) has been replaced by /. and in the subsequent the notation 
/ = /(to) will imply an element of J-jv- 

3.2. Permutation invariant form of equations 



The BBGKY hierarchy (|22p can be written in an explicitly permutation invariant 
form, which is more concise and facilitates further analysis. 

Define a cyclic permutation operator, acting on X = {x, y, z}, by 

n(x, y, z) = (y, z, x) (38) 

and three reflection operators r x , r y , r z which act by permuting pairs of elements. For 
example 

r x (x,y,z) = (x,z,y) (39) 

The group of permutations in three dimensions consists of six elements, g = 
{i,7r,7r 2 ,ri,r2,r3}, where i is the identity operator i(x, y, z) = (x,y,z). 
Define twelve vectors tf, labelled by i = 1, 2, 3, 4 and a 6 {x,y, z}, by 

tf = (-1,1,1), # = (1,-1,1), t\ = (1,1,-1). 

* v = (a. -[.()). ti = fn.ti.iV 

(40) 



if = (1,0,0), # = (0,1,0), if = (0,0, 1). 



if = (-1,0,2), t| = (2,-1,0), *g = (0,2,-1). 
if = (-1,2,0), # = (2,0,-1), tf = (0,-1,2). 

Note that t y = 7rff, t z = irt y = Tt 2 t x , and t% = r a t^. Using these definitions, the 
coefficient BBGKY hierarchy equations ([35]) and (|37ap - (|37b|) can be written as 



d T fn(m) = J2 K a {m na m n2a f n -i(m - i?) - TO a /n+i(m + f?)} 

+ ]T TU a (m™ - ^ 2a ){TO a /n-l(™ - t a 2 ) - f n+1 (m + ig)} 



ael 
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+ W a {m" 2a (mf a - l)/„_i(m - - m^m™ - l)/„_i(m - 

ael 

- m»[/ n+1 (m + t£) - /„ +1 (m + £)] } (41) 
3.3. Representation by a graph 

To study the properties of (|41[) it is convenient to use the language of the graph theory. 
A graph Q = (971, 21) is an abstract representation of two sets: The set 971 of objects 
called nodes m G 971 will be identified with the set of indices (f3"4"|) of elements of Fm- 
Certain pairs of nodes m and mi of the graph will be connected by links, called edges 
and denoted by (m, mi). The set of all the edges is denoted by 21. 

Couplings among tensor coefficients in (|4"T]) [or in (|24|) or (|35|) which are equivalent 
to (|4"Tj) ] will be allowed to determine 21 by the following definition. 

Definition 1. An edge (m, mi) is an ordered pair of elements m, mi G 97t. (m, mi) G 
21 if and only if two conditions are satisfied: d T f n (m) is coupled with / ni (mi) by |^i| ) 
and t/ie coupling coefficient is not equal to zero. 

Two more definitions will help establish the terminology. 

Definition 2. The node mi is said to be adjacent to m if (m,mi) G 21. A set 
adj(m) G 971 is a set of all of the nodes, adjacent to m. 

Definition 3. A value function ip : 21 — > R is equal to the coupling coefficient between 
drfnijn) and f ni (mi) and is denoted by tp[(m,mi)}. 

The three-term recurrence structure of the BBGKY hierarchy implies that 
adj(97t„) C {97t n _i, 97l n+ i}. The adjacency relations will be labelled by ± to 
distinguish this structure. 

Definition 4. A set adj ± (m) C 971 is a set of all the nodes mi such that mi is 
adjacent to m and \mi\ = \m\ ± 1. 

By this definition, 

adj ± (97T„) C97W (42) 

Having established the terminology of graph theory, let us now determine the 
implications of the structure of the BBGKY hierarchy (T4"Tj) on the structure of the 
graph Q. For this purpose the following two lemmas are provided. 

Assume that m = (m x ,m v ,m z ) G 971, mi = (mf , mf, mf) G 97t, mi G adj + (m). 

Lemma 1. Nodes m and mi are connected by a edge (m, mi) G 21 iff for some a £ I 
both conditions in one of J^ffap - f^gaj ) are satisfied 

mi =m + ti ip\(m,mi)] = —K a m a =/= 0. (43a) 

m 1 = m + t% ip[(m,mi)] = -W a (m™ -m^ a ) ^ 0. (436) 

mi = m + t% ip[(m,mi)] = -W a m a ^ 0. (43c) 

mi=m + i4 ip[(m,mi)] = W a m a ^ 0. (43d) 

Proof. By definition [T] and by assumption mi G adj + (m), an edge exists iff the term 
d T fn( m ) in dUl has a non-zero coupling with /„+i(mi). Since Va, i : \m±t?\ = 
\m\ ± 1, the first necessary condition for (m, mi) G 21 is that mi is equal to one of 
the terms listed in the first column of ()43a[) (|43<i|) . The second necessary condition is 
non-vanishing of the corresponding coupling coefficient, which are listed in the second 
column of (|4^o |) -([43 ^) . □ 
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Lemma 2. (m,TOi) G 21 implies (mi, to) £ 21. 

Proof. Assume that (to, TOi) S 2t. By lemma [TJ toi = to + 1" for some i = 1,2,3,4 
and some fixed a6l By definition^ the adjacent set adj - (mi) is determined by the 
non-zero couplings between d T fn+i{m\) and f n (m2) where to is substituted by mi in 
(|4ip . By (|4"Tj) and lemma [1] to 2 = m + <° — t*, therefore to 2 = to if i = j and a = b. 
Examining the coupling coefficients for each i and a in (|4"Tj) one finds 



V>[(toi,to)] = K a (m™ + l)(m? a + 1) i = 1. (44a) 

^[(m^m^^K + llK'-ra") i = 2. (446) 

^[(m 1 ,m)]=ft' (l (ffl A + l)(m A +2) i = 3. (44c) 

^[(mi,m)] = -F a K a + l)(m™ + 2) i = 4. (44 d) 



Denote the absolute value of an integer coefficient of ip[(m, mi)], multiplying a coupling 
constant K a or W a , by j/(to,toi), and the one of ^[(mijm)] by v(m\,m). By 
assumption that ip[(m, Toi)] 7^ 0, it follows that the relevant coupling constant K a ^ 
or W a ^ 0, and v(m,mi) ^ 0. Comparing (|43ap - (|43dj) with (|44ap — (|44<ij) and using 
to° + 1 > 0, it follows that under these conditions v(mi,m) ^ as well. □ 

Comparing (|43a|) - (|43dj) with ()44a|) - (j44dj) we conclude that v(m,m\) ^ v{m\,m) 
in general. In particular, 

_ i)[(m,m{j\ = v(m,m x ) g 
■0[(mi,m)] ^(toi,to) 
is equal to a ratio of two positive integers. For mi = m + 1", define 7? by 



7? = -^[(to, to + i?)]V[(m + i?, to)]. (46) 

The following formulae for 7? are derived from lemmas 
2 

7° = X a max (m x , m x ) max (to 27 , m^) max (m z , ml) , (47a) 

l2=W 2 a {m a + l)(m" a -m* 2 ^ 2 , (47b) 

7 £ = M^TO a (m 7r2a + l)(TO 7r2a + 2), (47 c) 

7| =W 2 a m a (m™ + 1) (to™ +2). (47d) 



The graph Q is drawn as points and arrows on a plane. Since edges are directional, 
an edge (to, mi), connecting the m-th and TOi-st nodes can be represented by an arrow, 
starting at the point to and ending at mi. However, by lemma [21 all arrows come in 
pairs. Therefore, instead of connecting each adjacent pair of nodes by a pair of arrows, 
each pair of arrows can be substituted by a single line and an integer number on each 
side of the line to distinguish the direction of a link. The following convention will be 
adopted. For a pair of adjacent nodes m and toi G adj + (TO,), v(m, Toi) is to the right 
of (or below), and ^(toi, to) is to the left of (or above) the line, connecting m and Toi. 

To take into account the adjacency relation (|42j) and the composite structure (f34"j) . 
all the nodes m G 9Jt„ with a fixed n will be drawn at the same height. For example, 
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a segment of Q 1 § is shown in (|48[) . where each node is denoted by m z 



.m y m 




5.02 7.00 



(48) 



1.00 



4. Initial value problem for the coefficient BBGKY hierarchy 



4-1. Generalized Liouville operator 

Time evolution of coefficients / is obtained by solving the initial value problem 

d T f(T,f i )=J?f(T,f i ), /(0,/i) = /i, (49) 

where Jz?/(t, /i) denotes the right-hand side of the BBGKY hierarchy (|35|) . The initial 
condition /j is determined from the initial condition for the operator fi.„jv- 

The linear operator Jz? : J 7 ^ ^ J~n 1S called the generalized Liouville operator [H] . 
It is represented by a real- valued matrix, which multiplies vectors / = (/i, . . . , /at) t . 
The latter consist of N blocks by the composition formula (|28|) . The n-th block of 
/ consists of all the coefficients in the n-spin vector space arranged in a sequence. 
This block structure of the vector space Tn imparts a similar block structure on Jzf , 
comprising N x N rectangular matrices L™ e ^ d (v xd (v j where 1 ^ i,j ^ iV. The 
three-term recurrence structure of (pt5j) implies that the only non-zero blocks L™' are 
such blocks that \i — j\ = 1. Define 



^(n+l.n) ^ ^d(n-\-l)xd(n) 



(50) 



The generalized Liouville operator Jz? has the form of a block tridiagonal matrix with 
rectangular off-diagonal blocks and L~ 

/ -L+ \ 



if 





L.7 



-i3 
o 



-L 



L 



N-l 



N-l 

o / 



(51) 



§ Some nodes and vertices of Q 1 are not drawn in order to keep the image clear. 
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4-2. Self-adjoining transformation 

The BBGKY hierarchy ([55)1 can be rewritten using the graph notation as 
d T fn(m)= ^2 ip[(m,mi)]f n -i(mi) + ^ ^[(m,mi)]/ n+ i(mi). (52) 

miGadj - (m) mi £adj+ (m) 

Comparing this equation with (|5"T|) we find that for any to 6 9Jt„ and mi <E 9Jt tl +i, 
one matrix element of is equal to —tf>[(m,mi)] and similarly, one element of L~ 
is equal to ip[(mi,m)]. Relations between the indices i and j of a particular element 
(L^)ij and value functions tp are determined by an arbitrary choice of a mapping of 
m e 9Jl„ onto a sequence of integers {1, . . . , d(n)} for each n. It is not difficult to see 
that (L^)ij = —ip\(m, mi)] for some i and j implies (L~)ji = ip\{rn\, to)] for the same 
pair i, j. By a corollary of lemma [2] it follows that in general (L+)ij ^ (L~)ji, i.e. 
that are not mutually adjoint, and as a result that J>f is not anti-Hermitian 

L~^L+ T -> JJf^-jSft. (53) 

The evolution operator is formally obtained by exponentiating Jzf , 

Sf(T) = e^ r . (54) 

As a consequence of (|5"3")l , S^(t) does not preserve the conventional inner product 
(/iff) = 2 /ifli °f /iff ^ -^v for infinitesimal r and therefore S^(t) is not unitary with 
respect to this inner product. To get a positive result, one could try a change of basis 
and a similarity transformation of „£? by an invertible D such that si ', defined by 

*f = D _1 JSfD (55) 

is anti-Hermitian, 

^ = -^. (56) 

Since for a tridiagonal matrix (a matrix of the form (|51[) where are scalars) this 
can be achieved by a scaling transformation, consider D, defined by 

D = diag(£)i, J D 2 ,...,D JV ), (57) 

where Di is a d(i) x d(i) square matrix which is yet to be found. The transformation 
(|55p by (|57f preserves the block tridiagonal structure of ££ . In analogy with (|50)) - ((5T|) . 
denote the off-diagonal blocks of si by H^- By substituting into (|55|) . the condition 
(|56p is equivalent to 

Hn = D n l L+D n+ i, H~ = D~ +1 L~D n . (58) 
The condition that jz/ is anti-Hermitian is 

H~ = H+ T (59) 
and from and (|59|) it follows that (fS"6"f is equivalent to existence of a solution of 

{D n+1 Dl +1 )L+ T = L-(D n Dl), n = l,...,7V-l, (60) 

where D n D^ must be determined by solving (|60[) starting from an arbitrarily chosen 
D\. For each n, ([6U)l is a set of <i(n) x d(n + 1) equations for the unknown elements 
of a square symmetric matrix, which has d(n + l)[d(n + 1) + l]/2 free parameters. 
A necessary condition for solvability of (|60j) is a greater number of unknowns than 
the number of equations, which amounts to satisfying [d(n + 1) + l]/2 ^ d(n). Using 
(|25p this inequality simplifies to 2n + 6 ^ which is certainly satisfied for n > 0. To 
construct an explicit form of si, the following theorem is provided. 
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Theorem 1. Equations i (55|) -i f57j)) are solvable by a diagonal matrix D. Each element 
of of srf is equal to ^J^if , defined by J^6| ) ; for some i and a. 



Proof. Assume c S J-pj such that c(m) ^ OVto G 9Jt. Consider a transformation of / 

f' n {m) = c{m)f n {m) (61) 

This transformation is of the form (|55p with a diagonal D, where = l/c(m) Vi and 
some m = m(i) e 9JT. It is convenient to express the BBGKY hierarchy equations for 
fn(m) in the form (|52|). Substitution of (|61j) into (|52p results in the equation of the 
same form as (|52|) with /s replaced by /'s, and ■0s replaced by ip's, and where the 
latter are defined by 

tp'[(m,mi)] =i)[(m,m 1 )}c(m)/c(m 1 ) (62) 

The proof will be established if c(m) can be determined from (|62"T) . and the condition 

tp'[(m,mi)] = -ip'[(mi,m)] (63) 

is satisfied for each to <E 071 and toi S adj(m), as this condition is equivalent to ([51)1) . 
Combining ([62]) and (|63|) . and using (|45|) one obtains 



Q . . O . . ^(TOl , TO) 

c 2 (to) = c 2 (toi ) - ( . 64) 
v(m, nil ) 

From (|64p it follows that c(to) is determined by its adjacent node mi. Equation ([54^1 
can be applied recursively to determine c(m) by c(mi), where toi is any node such 
that to is reachable from toi , multiplied by a ratio of products of numerical factors v 
along the path, connecting toi with to. 

Definition 5. A path p of length k > is an ordered sequence of k edges (mj,mj+i), 
1 ^ j ' ^ k such that rrij±i is adjacent to rrij, 

p(mi -> m k +i) = ((mi, to 2 ), (to 2 , m 3 ), . . . , (m k , m k +i)). (65) 

A loop po is a closed path p, i.e. a path where m^+i = mi 

Pa = ((toi, to 2 ), (m 2 , m 3 ), . . . , (m k , mi)). (66) 

Using this definition 

c 2 (m k+1 ) = c 2 (m l ) ; ' , , 67 

J -- L v(m" (e).m' (e)) 

where e denotes an element of p such that e = (m'(e), m"(e)). Since most to's have 
more than one adjacent to", most pairs of connected nodes can be linked by several 
distinct paths. Therefore, solvability of (|55"]) by a diagonal D rests upon the condition 
that c(TOfc-)-i) is uniquely determined by toi and is independent of p. Looking at a 
segment of Q 1 in (|48p it is clear that any path p(m\ — > mk+i) could be augmented 
by inserting loops at any of the intermediate nodes. Therefore, independence on the 
traversal is achieved if for any loop p 

n v{m'{e),m"{e)) = 
^f(ro"(e),ro"(e)) ' 1 ' 

This condition could be turned into an additive condition for loops by taking a 
logarithm of (|68p and defining the current j on each edge. 
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Definition 6. For any m! G 911 and m" S adj(m'), and e = (m', m") G 21, the current 
j(e) G E is defined by 

j(e) = log [i/(m', m")/v(m", m')] . (69) 

Solvability of ([55)1 by a scaling transformation can be formulated as a "zero 
circulation'' condition: For any loop po we require that 

E^ e ) = °- ( 7 °) 

eGpo 

Considering lines that make up any loop in (|48[1 as boundaries of that loop that 
separate the "inside" from the "outside", the zero circulation conditions is verified 
by multiplying all the numbers along the loop which are either outside or inside 
its boundary. As can be seen from (|4"5)l . the most elementary non-trivial loop po 
involves four nodes, and more complicated loops can be constructed from a sequence 
of elementary loops. It is therefore sufficient to verify the circulation condition ([70)1 
for each loop of the form 



m + t" ' m + t) 



where tf ^ tj. As an example, consider in (gSJ) a loop, that visits {102, 211, 122, 013}. 
We find n i Vi = 2 • 2 • 2 • 3 = 24, and n* = 2 • 4 • 3 ■ 1 = 24. Therefore £2) is satisfied 
for this particular loop. A straightforward calculation (not detailed here) confirms 
that ([71)) is satisfied for arbitrary tf ^tj. 

The elements of srf are expressed in terms of value functions tp'[(m,mi)]. 
Substituting ([6"4"| into ([52^1 . with k = 1, such that m,2 = mi + , and using (|46l) 

V>'[(m, m + i?)] = yj-ip[(m, m + t?)M(m + tf,m)] = ^7f . (72) 

It follows that ijj'[(m + tf, m)] = ip'[(m, m + tf)]. □ 

Each coefficient c(m) can be calculated by following the path to m from some 
fixed reachable mo, whose coefficient can be chosen arbitrarily. Setting c(mo) = 1, 

cH = ^»^{ n Tln m "n! } 1/2 ^ 0*) 

I - LJ - vim [e),m [e}} j 

where p(mo — > m) is any path, connecting mo with m. 

Since d T D~ 1 f = sz/D~ 1 f where is anti-Hermitian, and (D~ 1 f)i — c(m)f(m) 
for some m = m(j) by theorem [TJ we conclude that the evolution operator ([54)) is 
unitary, provided that the inner product on the vector space J-jv is defined as follows 

N 

</ 1 '/ 2 >^ = E E c 2 (m)^(m)/^(m) f\feF N . (74) 
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5. Invariant subspace decomposition 

The following definition of an induced graph will be useful for the discussion. 

Definition 7. A graph Q 1 = (9ft 1 , 21 1 ) = ind(</,mi) is a subgraph of Q = (271,31) 
induced by mi € 9JI, if 97T 1 and 21 1 are defined as follows. The set 9J1 1 C 271 consists 
of all the nodes m £ 971 such that mi is reachable from m, and 21 1 C 21 is the set of 
all the edges (m, mi)Vm, by which mi can be reached from m. 

A graph Q is fully connected if all of its nodes are reachable from any one of 
them. The same can be stated using induced graphs as ind((J, m) = Q for any m G 9Jt. 
Conversely, if for some m\ G 271, ind(C/,«ii) ^ G, then Q is a disconnected graph. A 
disconnected graph is a union of several connected graphs Q l : Q = Q 1 U Q 2 U . . . 

If the BBGKY hierarchy (f5S"|) is represented by a disconnected graph Q, then 
all expansion coefficients in J~n can be grouped by the node set 97T of its each 
connected subgraph Q 1 = (dJl l ,%l l ). Each group of coefficients evolves by a subset 
of BBGKY hierarchy equations for these coefficients. Since different subgraphs Q % are 
disconnected, the BBGKY hierarchy equations do not couple coefficients in different 
SETPs and therefore each group of coefficients evolves independently from the remaining 
coefficients. Thus each set of nodes 97P defines indices of a group of expansion 
coefficients that belong to an invariant subspace F l N of the full vector space Fn- 
In this section it will be shown that the connectivity of Q depends on the structure of 
the effective couplings, and therefore J determines the topology of Fn ■ 



5.1. Effective coupling constants 

In [6] and in section[2]it was shown that the effects of magnetic field h can be accounted 
for, in thermodynamic limit, by replacing the matrix of coupling constants J by a 
matrix of effective coupling constants J = (J a b) which is defined by (|23|). If h = 0, the 
interaction picture is not necessary since the equations with a non-averaged J are of 
the same form (|24[) . In this case assume J — J. A common, but not the most general 
setting is when the couplings J are diagonal, 

J = di&g( J x ,J y ,J z ) (75) 

and h = (h x , h v , h z ) is an arbitrary vector with a condition h ^= 0. Define h = h/ \h\. 
In this case, the matrix of effective coupling constants J can be written as 

where 

Jd = {hl(J y + J,) + h 2 y (J x + J z ) + h 2 z (J :c + J y )} l 3x3 (766) 

Js = [Jx + J y + Jz - 3(h 2 x J x + h 2 y J y + h 2 z J z )} {h a ® h b ) (76c) 

and where (h a <g> hb) denotes the outer product of h = (h x , h y , h z ). 

Inspecting (|76ap - (|76cp we conclude that if J is given by ([75)1 and either h = or 
h has only one non-zero element [for example h = (0, 0, h z )\, then J s = 0^x3 and J is a 
diagonal matrix. In general, if two components of h are non-zero, as in h = (h x , 0, h z ), 
then J has one non-zero off-diagonal element; if all three components of h are non-zero, 
then all off-diagonal elements of J are non-zero in general. The first two cases will be 
discussed in detail in the following, while it is assumed that Q is fully connected if all 
the diagonal elements of J are unequal, and all off-diagonal elements are non-zero. 
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5.2. Diagonal couplings 

In this section, a special case will be considered where the coupling matrix J is diagonal 

J = di&g(J x ,J y ,J z ) . (77) 

Define J to be anisotropic if all three coupling constants J a are distinct. Equivalently 
J is anisotropic if K x = J y — J z ^ 0, K y = J z — J x ^ 0, and K z = — J y ^ 0. A 
similar setup has been studied in [5] where the coupling matrix was diagonal but not 
anisotropic. This special case will be further studied in section [5^1 

Equation ([77)1 implies W x = W y = W z = 0. As a result, v^{m,f) simplify 
considerably 

v~(m; f) = K x m y m z f(m x + l,m y - 1, m z - 1) + K y m z m x f(m x - 1,7% + l,m z - 1) 
+ K z m x m y f(m x - l,m y - l,m z + 1) . (78a) 

v+(m;f) = - K x m x f(m x -l,m y + l,m z + 1) - K y m y f(m x + l,m y - l,m 2 + 1) 

- K z m z f{m x + l,m y + l,m 2 - 1) . (786) 

It follows from (f78a|) - (|78 6[) that all the edges of Q are generated by t\, defined by (|30]). 
Assume that J is anisotropic. Define 





{(1,0,0)}, 01? = {(0,1,0)}, 01? = {(0,0,1)}, 


(79 a) 




{(2,0,0), (0,2,0), (0,0,2)}. 


(796) 


;rapht 


; of induced by these sets 




^ = 


(gi 1 ,2i 1 ) = ind(e?,^), 


(80a) 


S 2 = 


(W 2 ,2l 2 ) =ind(^ J *n?) J 


(806) 


3 = 


(9T 3 ,2t 3 ) =md(S,0rtf), 


(80c) 


S 4 = 


(0t 4 ,2l 4 ) = ind(a,^). 


(80 d) 



Here m % n C £0T„. Similarly to the composition of DJP- N ' ([M]) . we define 



JV 



«n w) = |J mi . (8i) 



n=l 



Set mf is empty. The following lemma describes the structure of the remaining m i n . 
Lemma 3. For any 1 ^ n < N — 2, and i = 1, 2, 3, 4, 

<+2= (J {™+(2,0,0),m+(0,2,0),m+(0,0,2)} (82) 

Proof. Assume that n is fixed and that (|52"j) is valid for 71 — 2 and n. Induction of 
nodes in one step from m z n to m i n+l is obtained by following all the edges which are 
defined by vectors t": m — > m + if. Induction in two steps results in 

K+2= (J {"* + {*? + **}}■ (83) 

abel 

where 

{*? + A) = {(0, 0, 2), (0, 2, 0), (2, 0, 0), (-2, 2, 2), (2, -2, 2), (2, 2, -2)} (84) 

The first three terms in (|54"|) satisfy (|52"j) . Each of the remaining three terms can 
be viewed as a three step path, e.g. (-2, 2, 2) = (-2, 0, 0) + (0, 2, 0) + (0, 0, 2). By 
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assumption ml = m + (-2, 0, 0) G W n _ 2 and m" = m' + (0, 2, 0) G W n . The last step 
m'" = m" + (0, 0, 2) satisfies @ therefore m'" G «Jl* 1+2 . 

By direct calculation one verifies that 0^ = {(0,1,1)}, 9t£ = {(1, 0, 2), (1, 2, 0)}, 
and = {(2, 1, 1), (0, 1, 3), (0, 3, 1)} satisfy the assumption of the lemma for n = 1,3, 
and n = 2,4. The corresponding sets with i = 2, 3 are obtained by permutations of 

The case of i — 4 also satisfies the assumption with n = 2, 4 and n = 3, 5 as 9t| = 
{(1, 1, 1)}, mi = {(0, 2, 2), (2, 0,2), (2, 2,0)}, and (ttf = {(3, 1,1), (1,3,1), (1,1, 3)}. □ 

The main result of this section is contained in the following theorem. 

Theorem 2. Assume that J is a diagonal anisotropic matrix. Then Q % are disjoint 
subgraphs ofQ, and Q = G 1 U G 2 U G 3 U G 4 . 

Proof. For the first part of the theorem it is sufficient to show that for each 1 ^ n ^ N, 
9^n9T£ = i^j,Vn. (85) 
For the second part one must show that 

4 

Win = |J K Vn. (86) 

i=i 

In analogy with (|25p , denote the number of elements in by di (n) . The connectivity 
structure of Q 1 , G 2 and G 3 is the same because their node sets W are related by cyclic 
permutations, in addition to K x ^ 0, K y ^ 0, K z ^ 0, implied by the anisotropy of 
J. Therefore the corresponding node sets 91^ have the same number of elements for 
each n, and the condition (|56")) can be replaced by (|55|) together with 

d(ri) = 3di(n) + d 4 (n) Vn. (87) 

For n = 1,2 conditions (|85[) and (|86| are verified directly. In particular 

OTi = {(1,0, 0), (0, 1, 0), (0, 0, 1)} = uiiWr , (88a) 

M 2 = {(2, 0, 0), (0, 2, 0), (1, 1, 0), (1, 0, 1), (0, 1, 1)} = Uii^ . (88b) 

Consider the parity of m x , m v , m z of the elements m = (m x , m v , m z ) G 9^. Induction 
of nodes in one step by m — > m + t", adds or subtracts unity to each element m° of m. 
This operation inverts the parity of each m a in one step. In two steps, induction does 
not alter parity, as can be seen from lemma [3] By examining OT^, 91 2 and applying 
this simple parity alternation rule, the parity composition of 0^ can be determined. 
Define i a by i x = 1, i y = 2, i z = 3. For all k = 0, 1, . . . , N/2 and a = x,y,z 



m G yC 2 a k : m a is even, m Ka , m 7r a are both odd, (89a) 

m G 9T 2 fc+i : TO ° is °dd, m7rQ ' m ^ a are both even - (^96) 
Similarly, if i = 4, 

m G 0T^ fc : m 21 , m y , m z are all even, (89c) 

m G ^2fe+i : mX i m% ') mZ are au °dd. (89d) 



It follows from (|89a[) — i|89<ij) that the patterns of parity are incompatible among all Vl^ 
for each n, and therefore these sets are disjoint. This proves (|55")) . 



Symmetries of a mean-field spin model 



18 



Counting of the elements in 71^ and can be achieved by parameterizing all 
m = (m x , m v , m z ) with the necessary properties (|89ap - (|89rf}) . Define 

m kl3 ; = (2k-2i,2i-2j,2j) k^i^j^O. (90) 

Using this definition all the cases in (|89o,[) (|89(i|) can be parameterized as follows, 

^k+i = {m = mkij + (1,0,0)} fc^O, (91a) 
^2fc+2 = {m = m fc0 - + (0,1,1)} fc^O, (916) 

^2fe = { m = m kij} k ^ 1 , (91c) 

^2fe+3 = {m = m WJ + (1,1, 1)} fc^O. (91 d) 

Parametrizations of *Tt^, 9t 3 are derived similarly from 91*. From (|91ap - (|91<ij) it follows 
that the number elements in various sets Vl l n is related to the number of elements myt 
with a fixed k by 

di(2Jfe + 1) = di(2fe + 2) = d 4 (2fc) = d 4 (2fc + 3) 

it 

= E C 1 **^) (Wo) = E( j + !) = ■ ( 92 ) 

i,3 i=0 

where l_y = 1 if A" is true and lx = if X is false. Counting of terms by (|91a[) - (l91«D . 
done separately for the even n = 2k and odd n = 2k + 1 results in, respectively 

3di(2fc) + d 4 (2k) = M(k - 1) + d(k) = d(2k) , (93a) 
3di(2fc + 1) + d 4 (2/c + 1) = 3d(fc) + d(fc - 1) = d(2fc + 1) . (936) 

Since the argument on the right hand side is n in both cases, this proves (|87[) . □ 
5.5. Non-diagonal couplings 

Assume that one off-diagonal element of the coupling matrix is non-zero, i.e. W a + 1 
for some a £ {x, y, z}, and W na = W 7T 2 a = 0. For example if a = x, J has the form 



J X 


w x 





w x 













X 



Define 

y\ A = W a urn 4 , m B = u of* 2 « . (95) 

Define graphs t?' 4 C (? and Q B C 5 as induced from these sets 

A = (9t A , 2l A ) = ind(£, m A ) , (96a) 
B = (9t A ,2l A ) = ind(g,9? B ). (966) 

Theorem 3. Assume that K x ^ 0, K y ^ 0, A' z ^ 0, and W a ^ for some 
a E {x, y, z} but W^a = W n 2 a = 0. Then graphs Q A and Q B are disjoint subgraphs of 

g, andg = g A ug B . 

Proof. Here, the complete node sets Ot" 4 and 91 s are provided and it is sufficient 
to show that induction of nodes, generated by vectors tf, t\, if, and t%, if, t% 
which correspond to non-zero effective coupling coefficients, leaves the node sets 
invariant. The first three vectors generate the sets Ot 1 / 2 / 3 / 4 , discussed extensively 
in section 15.21 From the definition P0|) it follows that induction in one step by 
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the remaining three vectors inverts the parity of m a only, and does not change 
the parity of m™ and m 71 " a . By inspecting the parity structure of 91^, discussed 
in theorem [21 it follows that the addition of these vectors modifies the adjacency 
relations to adj ± (W„«) = {9^, 9t* ±1 }, adj ± (^ i ) = {9t* ±1 , 9^}, and similarly, 

adj ± ({D r t* l ,ra ,9 r tJi ,r2 °}) = {91^!, from which the invariance requirement of 
Ot" 4 and 91 s follows. Furthermore, disjointness and completeness of the graphs Q A 
and Q B follows from from the corresponding properties of Q l and ®iQ l . proved in 
theorem [5J □ 

5.4- Isotropic diagonal couplings 

A rather different situation is encountered when J is diagonal and isotropic, as in 

J = diag(J ± ,7 ± , J z ). (97) 

Now we have W x = W v = W z = K z = 0. Since ^ - Ki = 0, K must be of the form 

K = (-k,k,0), (98) 

where k — J, — Jj_. Since n is the only coupling coefficient appearing in each term 
on the right-hand side of (|4T]) . it can be scaled to ±1 by an additional transformation 
of the time variable. Since the evolution is unitary, we may set k = +1, and thus 

d T f n (m) = m x f n+1 (m x - l,m v + l,m z + 1) - m y f n+1 {m x + l,m v - l,m z + 1) 
+m z [m x f n - 1 (m x - l,m v + l,m z - 1) - m v f n -i(m x + l,m y - l,m z - 1)] . (99) 

Since each term on the right-hand side of (|M]) is multiplied by m x or m y , the left-hand 
size of (|99p is zero whenever m is of the form m„ = (0, 0, n). As a result, f n {m n ) is a 
constant of motion for each n = 1, 2, . . . , N. 

To determine the remaining invariant subspaces, it is instructive to study a graph 
that corresponds to d%f = -Sf 2 /. A straightforward differentiation of (|9"9"]l results in 

d 2 T f n {m) = - (2m 2 + l){C f„(m) - C x f„(m + t x ) - C y f n (m - t x )} 

+ C x {f n+2 (m + t ± +t ll ) + C z f n - 2 (m + t ± -t ll )} 

+ C y {f n+2 (m -tj_ + t\\) + C z f n - 2 (m -t±- t\\)} 

~ C {f n+2 (m + + C z f n - 2 (m - t||)} . (100) 

where 

C a {m) = m a (m a - 1) ae{x,y,z}, C (m) = 2m x m y + m x + m y , (101a) 
and where 

ii_ = (-2,2,0) i|, = (0,0,2). (1016) 

As in previous sections, we associate a graph Q = (9Jt, 21) to the operator J§f 2 on the 
right-hand side of (|100|) . and determine the disconnected subgraphs of Q by following 
all possible paths, induced by f|101 £>[) from various initial elements m. For example, all 
changes in m z are provided by ±tj_ , which implies that the node set of any disconnected 
subgraph Q l C Q is such that m z is either an even or an odd integer from an interval 
[mfnin, N], where ^ "m-min < N is the smallest possible value of m z of a particular 
subgraph. Moreover, if Q l is a disconnected subgraph with m z = 0, 2, . . ., then another 
graph Q n , which differs from Q l only by replacing m z by m z + 1, is also a disconnected 
subgraph of Q: Q h C Q. 
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The "horizontal" induction of nodes (by paths that do not change m z , nor 
\m\) is obtained by applying m i— > m ± t±_ iteratively with the usual constraint 
m G 93t. Four cases are possible, depending on the parity of m x and m y . Define 

Ol,kn, 02,kn Q 9Jl„+2Jb+l, 03,fc« !_ 9Kn+2fc- an d 04,fcn C 9JT„+2A:+3 by 
fc fc 

oife, = (J (1 + 2j, 2fc - 2j, n) , o 2fcn = U (2* - 2j, 1 + 2j, n) , (102a) 

3=0 j=0 
k k 

o 3 fc» = |J (2k - 2j, 2j,n) , o 4fc „ = (J (1 + 2k - 2j, 1 + 2j,n) . (1026) 

3=0 j=0 

Each of the sets 0ifc„ C DJt n i is invariant with respect to the horizontal induction. 
Define even and odd node sets by combining the corresponding 0jfc„ as follows 

N/2 N/2 

O e lk =\J o lk . 2n D° k = [J ol >2n+1 . (103) 

n=0 n=Q 

By construction, O^l" are invariant with respect to node induction by (|1016[) , therefore 
a family of induced graphs, parametrized by e/o, i = 1, 2, 3, 4 and fc = 0, 1, . . . , iV/2 
and defined by 

<3t = (StXl°) = ^(g i Dt) k = 0,1,..., N/2, i = 1,2,3,4 (104) 

are disconnected subgraphs of (J. 

The range of the parameter k in (| 104)) implies that the number of disconnected 
subgraphs of £7 scales as O(N) and consequently, the same scaling holds for the number 
of invariant vector spaces J~^l° N , associated to each disconnected subgraph Q\l° ■ 

Associated to each 0i kn there is an n'-spin vector space S n i , where n' = n + 2k + 1 
if i = 1,2, ri = n+2kiii = 3, and n' = n.+2fc+3ifi = 4. By definition ([Waj) - (fT02b|) . 
Oikn has fc + 1 elements, thus each S n i subspace is fc + 1 dimensional. As a result, the 
operator _S? 2 associated to a disconnected subgraph Q^l° has a block structure where 
all the blocks are square (fc + 1) x (fc + 1) size matrices, and where fc is fixed. 



5.5. Dependence on initial conditions 

In the preceding sections it was demonstrated that, depending on the structure of J, 
connectivity properties of Q vary. They may be summarized as follows. 

Q = |J G a * , (105) 

where a = {a>i, a<i, . . .} is a set of partition indices and where Q ai = (9J? Qi ,2l Qi ) is 
a connected graph for each on G a. If J is a diagonal anisotropic matrix as defined 
in section T5.21 a = {1,2,3,4}; if J is diagonal anisotropic and, in addition, one off- 
diagonal element of J is non-zero, as defined in section [S~3l a = {A, B}. Lastly, if J 
is a diagonal isotropic matrix then, as discussed in section IS~4l a — U k Uj (e/o,j, fc). 

A collection of components f(m), m G 9Jt Qi was viewed as a vector in a vector 
space J 7 ^ C /F N . Connectendess of a graph Q ai implies that the subset of BBGKY 
equations for the evolution of vectors in F 1 ^/ involves only vectors in J 7 ^ and therefore 
J 7 ^/ is an invariant subspace of J-jy. It follows that 

se = . (106) 

q 2 £a 
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where Jz? is the generalized Liouville operator such that Jzf Qi : n> J 7 ^/ acts on 
each invariant subspace. 

Almost every Q ai (the only exception is a — (e/o, 3, k) in section !?^)) slices across 
the entire graph Q, meaning that Q ai describes tensor coefficients from n-spin vector 
subspace with all n m ; n ^ n ^ N (possibly n is constrained to either even or odd 
values) where n m - m is some smallest order spin vector space <Sn min of that subspace 
•Fjv • As a result, the initial value problem (j49|) in each J 7 ^ has to be solved to all 
orders in the expansion coefficients f n in general. 

From (|106[) it follows that dynamics a-priori depends on the initial conditions. 
For example, if the initial condition is such that the only non-zero components of / 
are in one of the invariant subspaces J 7 ^ , then the properties of equilibration and the 
equilibrium state will be determined by, respectively, the spectrum and the null-space 
of _Sf Qi . Both these properties may depend on and thus on the initial state. 

6. Conclusions 

Aspects of the dynamics of the anisotropic Curie- Weiss quantum Heisenberg model 
are addressed by studying the structure of its equations of motion. This study is 
based on the BBGKY theory and a special expansion of reduced density operators in 
a basis of Pauli operators, and it builds on a recent work where the coefficient BBGKY 
hierarchy was derived [BJ. 

This work was motivated by a special case of coupling constants [BJ , where time 
dependence of a certain subset of spin-spin correlators could be determined analytically 
to all orders, allowing to determine the evolution of a class of observables. It was hoped 
that the method of solution could be extended to all observables, initial conditions, 
and more general coupling constants. As in became apparent that a special symmetry 
of the coefficient vector space was behind the analytic solution, a comprehensive study 
of symmetries of a general case was called for. 

To examine the structure of the coefficient vector space, the BBGKY hierarchy 
was represented by a graph, whose vertices are the expansion coefficients, and edges 
are couplings in the BBGKY hierarchy. It was shown that, under certain conditions, 
this graph is disconnected. A disconnected graph implies that certain groups of 
expansion coefficients are causally unrelated. To each such group one can associate an 
invariant subspace of the underlying vector space. It therefore follows that the vector 
space is partitioned into invariant subspaces, and that this partition depends on the 
structure of the matrix of effective coupling constants. Three scenarios of partition 
were identified in addition to the most general case: If J is a diagonal anisotropic 
matrix, there are four invariant subspaces. If in addition J has one non-zero off- 
diagonal element, there are two invariant subspaces. In the most general case, the 
vector space is irreducible, which corresponds to J that is diagonally anisotropic and 
has two or three non-zero off-diagonal elements. Lastly, it was found that, if J is 
diagonal and isotropic, the number of invariant subspaces scales as O(N) with the 
size of the system TV. 

In addition it was shown that, to make the evolution in the coefficient vector space 
unitary, it is necessary to scale the components of vectors by certain weight factors, 
which is equivalent to redefining the conventional inner product, where the products 
of vector components are multiplied by the squares of these factors. 

This work provides several insights about symmetries of the Hilbert space of a 
simple quantum model. It is believed in this work that similar symmetries are present 
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in a class of more realistic lattice spin models, where interactions among particles 
depend on the distance between them. It is also believed that these symmetries 
should be taken into consideration when low-dimensional approximations, such as 
kinetic theories, are constructed. 

Appendix A. Derivation of (|20a|) and (|20fe|) 

Using permutation symmetry of the tensors f%, components of defined by (|20ap - 
(|20fr[) can be written as 

n 

Kt) (ai - a " ) = E E £a<w (- J " w )/«+r 

uvw^X i—1 
n 

(„-)(«!...«») = E E e a > uv (-J a nfn-\ (Lit) 

where the dots on the right-hand side denote the multiindex (ai . . . a n ) modified by: 
deleting a, in (jl.la[) . and deleting both cij and aj in (|20a[) . 

For a S {x, y, z}, define b a = ira and c a = it 2 a. Define the permutation operator 
7r Q acting on a triple (a, o a ,c a ) by 7r a (a, 6 a ,c a ) = (x,y,z). Due to the permutation 
invariance of f n , the sums over i and j in (JTTTaJ) — (jTTTfoJ) can be replaced by sums 
over coordinate labels, multiplied by the number of occurrences of that label in 
the multiindex (ai...a„). For example in the term (| 1 . 1 a|) this multiplicity factor 
is equal to m a , by definition (|3"Tj) . Passing to a notation fn — > fn(m) in the term 
(t>+)( ai >"'' a ™) and performing u and w summation explicitly, this term is rewritten as 

(«+)(«*•••->= e ™^ a ™(-^)/« + Y 

= E m a {(-J baW )f n+ i[n a (m a -l,m ba +6(b a ,w),m Ca +1 + 5{c a ,w))) 

aw£T 

+Jc a wfn+i[Tr a (m a - 1 + 6{a, w) , m b " + l + 5(b a ,w),m c + S(c a ,w))]} (1.2) 

where 5(a, b) is the Kroenecker symbol, acting on coordinate labels. In arriving at this 
expression, note that the term m a — 1 arises because a,; is deleted from the expression 
/ra-i") while the terms m ba + 1 and m Ca + 1 arise from the w-summation. 

Performing the iu-summation, and passing to notation (i>+)( a i--- a ») — > ii+(m), we 

find 

v+(m) = E™ a { - (J baba - J CaCa )fn+Mm a - l,m b " + l,m c " + 1)] 

a 

-Jb a c a {fn+iK(m a - l,m\m c ° +2)] - /„ +1 [i fl (m" - 1, m K + 2, m c °)]) 

-4c„K» -m c =)/„ +1 [7r Q (m a + l,m 6 ",m c ")]} (1.3) 

To arrive at this expression from (|1.2p . a permutation invariance property was used: 
all the labels a, b a , c a can be cyclically permuted without changing the value of the 
sum, and moreover f[ir a (a, b a , c a )] = f[^b a (b a , c a , a)]. For example 

^2m a J Caa f n+1 [ir a (m a ,m ba + l,m c °)] = E m ° a Jb a cjn+i[^a(m a + l,m ba ,m Ca )] . 
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The final expression (|37a[) is obtained from (|1.3[) by identifying 

K a = J baba - J CaCa W a = J baCa (1.4) 

The term (|l.lb[) is treated similarly, i.e. the sums over i and j are replaced by sums 
over u, v € I. Since the labels and a,j may be equal for i ^ j there are two cases 
of multiplicity factors to consider, m ai (m ai — 1) if a, = aj, and m ai m a: > if a; 7^ a,. 
Separating these two cases, v~ is written as 

K7) (ai - a " ) = ]T rn a (m a -l)e auv (-J ua )fZz\+ ^ m a m w e auv (- J uw )f^\ 

auveZ auvw 

= ]Tm a (m a - l){4a/n-iKK - 2,m b » + l,m c °)] 

~J ba afn-ilKa(m a -2,m b \m c « +1)]} 

+ ]T mV{j Co J„-iKK - I,™ 6 " + 1 - %,»),™ c ° " «(ca,w))l 
aex 

w£{b a ,c a } 

~Jb a wfn-i[^a{m a - l,m b " - 8(b a ,w),m Ca + 1 - 8(c a ,w))}\ (1.5) 

In the second sum, performing the w-sum and after relabelling terms in the remaining 
a-sum, one finds after passing to notation {v~) a ^ n ^ — ¥ v~(m) 

v-(m) = J2{( jbaba - J CaCa )™ ba ™ Ca fn-il7Ta(m a + 1, m K - l,m c « - 1)] 

+m a (m a - l){JacJn-Mm a - 2,m K + l,m c «)] - J aba f n -Mm a - 2,m\m c " + 1 
+J baCa m a (m b ° -m c ») f n - 1 [7T a (m a - l,m\m c «)]} (1 
The final form (|376[) is obtained from (|1.6|) by substituting (|1.4[) . 
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